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Abstract 

It is shown that in presence of URCA shell pulsational energy losses due to neutrino emission 
and nonequilibrium beta heating are much less than energy losses by excitation of short- 
wavelength acoustic waves. Convective motion in presence of URCA shell is considered, and 
equations generalizing the mean free path model of convection are derived. 

2 Introduction 



It was suggested in |]T5[ that in the convective region cooling of matter may be enhanced in 
presence URCA shells appearing when matter contains an isotope with a threshold Fermi 
energy for an electron capture, coresponding to a density less then the central one. Presence 
of such isotope leads to existence of a jump in the composition at a density, corresponding 
to a threshold energy. During convective motion the matter in eddies around this density 
crosses periodically the boundary. That implies continious beta capture and beta decay in 
the matter of these eddies. 

Because of heating of a degenerate matter due to nonequilibrium beta processes , with 
account of convective URCA shell different conclusions had been done with respect to stabi- 
lizing or destabilizing the carbon burning in the convective degenerate core |8|, |l6|, ^, |l^, |ll|, 
12, H, 0, [l^. Here we calculate damping of stellar oscillations in presence of URCA shell (see 
analyze physical processes in the convective URCA shells and formulate approximate 
quantitative approach to the solution of this problem. 



3 Linear oscillations of a slab with a phase transition 

Consider the plane-parallel layer in the constant gravitational field with an acceleration g, with 
a phase transition at the pressure P^,, and polytropic equation of state P = Kp'^ , with K = Ki 
at P < and K = K2 at P > P,, Ki > K2. In static equilibrium P = gM {I - ^) , p = 

(^) ~ Tl)^^'^ ■ Here the pressure is continuous, but the density p has a jump at due 
to the jump of the constant K, M is the mass of one cm^ of the slab, m is the mass of one cm^ 
of the slab under the layer with a Lagrangian coordinate x, which is related to the density as 

p = (C — gx)^i-^ . In presence of a phase transition we have pQ = ■, Po = 

gM, gx, = C2-^-^, ^o = f, C2 = ^igM)'-^Kl/\ a = C2 + ^^{X - I). 
Here po and Pq are the density and the pressure at the bottom of the slab, xq and are 
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total thickness and thickness of the inner denser phase layer of the slab, A = p2*/pi*- The 
phase transition in the slab happens only if its specific mass M > P^/g = M^. 

Linear oscillations of the slab are reduced to Bessel equation for perturbations P,v = 
Pa,Vaexp{—iu!t) with solutions 



Pa = AV^J^ir]) + B^Y^irj), 
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where for two phases t]2 = ^^IoPq ' ~ ^2VA, = 1 
oscillations uj and relations between constants Ai, Bi, A2, B2 are obtained from boundary 
conditions and from relations on the phase jump ^ |l^]. The dispersion equation are 
obtained analytically for a frozen 
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and equilibrium phase transition 
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In the limiting case = 0, when the boundary between phases is on the inner boundary of 
the layer, and z^, = 1 the dispersion equation is reduced to J 1 (r2\/A) = for a frozen, and 



V\J^{n^/x) - (A - 1)^^ — -nj^{ny/x) = o 



(5) 
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for an equilibrium phase transitions. At m^: M, z^, ^ when the level between phases is 
moving to the outer boundary, we have the dispesion equation J 1 (0) = in both cases. 

Here 
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n = ^^, (6) 

To investigate the dependence of oscillattion modes with an ideal phase transition on 7 it is 
convenient to introduce i). = 2^^' with the first root — > at 7 ^ 00. All other roots 
tend to infinity at 7 ^ 00. First three roots of equations (^), are presented in Figs.l,2|^. 

Caption to Figures 1,2,3 

Frequencies Q, of the slab oscillations as functions of of the basic mode (Fig.l), modes 
with one (Fig. 2) and two (Fig. 3) nodes. Upper and lower curves correspond to frozen and 
equilibrium cases, relatively. Frequencies of these two cases coincide, when one of the node 
coincides with the phase transition. 

^Numerical solution of dispersion equations had been done by O.V.Shorokhov. 



3 



4 Damping of oscillations due to URCA shell in a highly de- 
generate matter 

Consider ultrarelativistic degenerate electron gas a good approximation in most URCA shells. 

n / 2 \4/3 

It corresponds to the polytrope with 7 = 4/3, and constant K = -^i^ sjlim^) ' '^here 

[^z = {^z\ + xz-i^^^ is the average number of nucleons on one electron. Here a two 
component mixture is considered consisting of elements with an atomic weight A and atomic 
numbers Z and Z — 1, with a beta transitions between them, xz and xz^i are mass concen- 
trations of these elements, xz + xz-i = 1. 

Let Ufe and 5 be Fermi energy plus rest mass energy of the electrons, and threshold energy 
for a beta capture, in units of mgC^; gz and gz-i be statistical weights of the elements {A, Z) 
and {A, Z — 1); -Fti/2 be a nondimensional value measured in the beta-decay experiments, or 
estimated theoretically. For small difference |(5 — nje| <C 1 we have simple expressions for 
an entropy increase during beta decay and capture in a fully degenerate matter 



(7) 

Here the rate of the entropy increase is equal to 1/3 of the energy loss rate by neutrino 
emission. During linear oscillations the beta reactions take place only in a thin layer of 
matter, crossing in its motion the boundary x = x*, m = m*, P = P*. In this layer the 
pressure may be represented by an expansion 

P = P^- g{m-m^)+P, P^'^ = P}'^ + ^P~^'\P - g{m - m^)], (8) 

Equations describing change of concentrations during oscillations averaged over the layer 
Am = Pa/ g, are written as 

dt Zbb gz 

at Ibb 

Here R = In 2g^;j^y^y^-^ ( ^^l^ ) • derive an equation describing decrease of the pulsa- 
tional amplitude, we should take into account the change of pressure due to change of the elec- 
tron concentration. Let p be a perturbation of the pressure determined by the beta reactions 
P = 5(M-m) + P + p, withp= Peg<P„P>P,; p = -IP^, P,g> 

P^ , P ^ -P* . 

In presence of damping we represent the velocity perturbation in a form v = Vait) sinwt, 
and define V{t) = Va out{'rn*,t)- so that v = V sintot. The amplitude of the perturbed outside 

'V^ 

P„ cos LOt 



pressure my be written as a function of y as Pa = zi ^^^p° /^f!!^*! ^; 7 = 1- The equation 



of motion gives the relation describing damping of oscillations in the layer Am(t) ^ , 
after subtraction of the proper oscillations of the slab at the frozen composition and averaging 
over the motion of the whole slab 
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dV An Am gp(m^) Am 

smut—— = — - — TT^r^r ~ itt^ (10) 

dt Am{t) M cos ujt M ' 



where Am = Taking into account, M = gPo, we get 



. dV gp{m^) , . 

sinujt—— = — . (11) 

dt Fq cos LOt 

Averaging over the oscihation period we obtain equation for decreasing of the amphtude of 
oscihations in the form 

'J=-Dv^ n=^(,,i^)i,xp,^r-^jp^^. ,12) 

dt 768ujPoZ V gz J JUvu) 

5 Energy balance and damping of oscillations 

Averaging equations (^) over the time and space we get an expression for heating rate of the 
oscihating slab due to nonequilibrium beta processes in the layer around the URCA shell 

Q^.U.^^^J,,SI^). R.^^.2 jf-f''^ 1^. (13) 

g 32x7571 \ gz J l2{Fti/2)z-iAmp m%c^ 

The rate of the neutrino energy losses Ly (ergs/s/cm^) is obtained by averaging over a time 
of losses in the whole oscillating layer. We get similar to ( p!^ ) 

Ly = 3Ri ^ Rsfl + ^V (14) 

g32 x 757r ^ \ gz ) ^ ^ 

Similar dependence of neutrino energy losses during oscillations because of URCA shell had 
been obtained by Tsuruta and Cameron (1970), who also took ( [T^ ) for the rate of loss of 
kinetic energy of oscillations. In the approximation of strong degeneracy the matter will be 
heated during oscillations with the rate (|l^), and neutrino luminosity is determined by (14). 



The source of both kind of energy fluxes is the pulsation energy of the slab, giving the rate 
of pulsation energy losses directly connected with beta reactions as 

MP) - in , r ^ _ ^^P*^ p5 (i , 9Z-1 

600^5 " V gz 



= -(Q. + L,) = -^^P'a ( 1 + ^ ) . (15) 



Defining E'pui ~ ^MV"^, and using approximately V = "^/^^j for URCA shell in the middle 
of the slab, we get from ([l^ ) the equation for decreasing of pulsational energy, connected with 
hydrodynamic processes, in the form 

pp-5/4 

p(dyn) _ p4/-| I gZ-1 . 

A ratio of damping rates of oscillation on i-th mode is . ^g"^) ~ ^ ^7 ■ follows that the main 

pul 

source of damping of oscillations in presence of URCA shell is connected not with the neutrino 
emission / nonequilibrium heating, but with a dynamical action of the nonequilibrium layer 
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of the slab, where beta reactions take place. This action leads to an exitation of short- 
wavelength acoustic waves with length / ~ <^ xq. When the wavelength of the excited 
eigen-mode approaches the thickness of the nonequilibrium layer Ax, formed by oscillations, 
both mechanisms of damping become comparable. For li ~ Ax = with account of relations 

^1 ^ ^0 ~ ^» ^ JV?' ^'^^ P* ~ Po, we get eI^JE^pJi''^ ~ 1. Importance of 

the dynamical damping of stellar oscillations in presence of URCA shell is connected with 
non-linearity of weak interaction rates, and deviations from eigen-oscillations under an action 
of weak interactions in this layer, leading to excitation of acoustic waves. The dissipation 
of pressure oscillations, connected with excitation of sound waves is inherent to any kind 
of dissipation, when the main term is nonlinear. The linear mechanism, connected with a 
conventional bulk viscosity, does not change the form of the eigenfunction of oscillations, and 
so no additional waves are excited. 



6 URCA shell in a convective motion 

It was concluded in Q that nonequilibrium heating is balanced by the change in convective 
flow, leading to the net cooling due to convective URCA shell. Nine years later same authors 



1 17 1 changed their mind, concluding that "convective URCA process can reduce the rate of 
heating by nuclear reactions but cannot result in a net decrease in entropy, and hence in 
temperature, for a constant or increasing density." This conclusion, as well as opposite one, 
made by using thermodynamic relations only, seems to be not convincing. Following this line, 
let us present two plausible scenarios, leading to two opposite conclusions. 

A. Due to action of a nonlinear bulk viscosity, the convection is damping in the vicinity 
of the URCA shell, decreasing the convective heat flux from the central part of the star. In 
the general heat balance of the star it means that cooling become less effective, and nuclear 
reactions become thermally unstable and lead to a nuclear explosion earlier, than without a 
presence of the URCA shell. Nonequibrium heating give additional heating, supporting the 
earlier nuclear explosion. 

B. Due to action of a nonlinear bulk viscosity, the convection is damping in the vicinity 
of the URCA shell, decreasing the convective heat flux from the central part of the star. Due 
to local decrease of the heat flux from the core the average temperature gradient increases, 
leading finally to increase of the convective flux soon after entering an URCA shell into a 
convective zone. If the increase of the convective flux prevails the nonequilibrium heating in 
the URCA shell, the general heat balance would be shifted to a larger temperature with more 
effective cooling, and the boundary of the thermal explosion would be postponed in time, if 
not eliminated. 

I cannot choose between these two scenario without construction of the numerical model 
taking into account all processes mentioned above. In such a highly nonlinear system, as 
a star with nuclear reactions, neutrino losses, degeneracy, convection, and many feedback 
influences it seems to be impossible to make a conclusion about the direction of process under 
the action of addtional URCA shell, basing only on thermodynamical ground. 

Convective modes belong to g-mode family, in which the local pressure perturbations are 
small, and could be neglected when the convective velocity is much less than the velocity of 
the sound. In this situation the sound wave dissipation of the convective modes imposed by 
URCA shell is negligible. The equations of stellar evolution in presence of URCA shell should 
take into account the following physical processes 
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1. loss of enegry due to neutrino emission in the URCA shell 

2. heating of the matter in the convective region around URCA shell due to nonequilibrium 
beta processes. 

3. Decrease of the convective velocity in the layer around the URCA shell due to energy 
dissipation connected with the nonequilibrium beta processes. Kinetic energy of the convec- 
tion is the source of energy for neutrino losses and nonequilibrium heating of the matter. In 
the condition of static equilibrium only energy and heat transfer equations should be modified. 
In the energy equation 

d5_(lE_P(fp_ _ , cu 1 dLr 

dt~ dt p^dt'^" '"^^'^ A^pr^dr' ^ 

in addition to other neutrino cooling processes e^, the new term is connected with heat- 
ing due to nonequilibrium beta processes around the URCA shell, having in mind strong 
degeneracy of electrons in this region. Neutrino emission in nonequilibrium URCA processes 
is accompanied by heating at high degeneracy, because the positive term pidn exceed the 
energy carried away by neutrino (Bisnovatyi-Kogan and Seidov, 1970). Convective motion, 
consisting of convective vortexes around an URCA shell is a source of additional neutrino 
energy losses, and of heating of the matter. This dissipation of convective energy may be 
described in the same way, as corresponding dissipation and heating during stellar pulsations. 
Therefore we use for description of these processes the formulae from the previous sections. 
If we accept that the pressure difference in the convective vortex is about one half of the local 
pressure, we use for the amplitude of pressure pulsations in (13) and (^) Pa = ^-P*. Taking 
also approximately Pq = -P*) Po = P*i and Uf^ = 6 we get 



Q^ = R^ , P = ln2^^^^ '—^ ^ 1 + ^^^ . (18) 



Equation (|1^) is related to energy losses averaged over the whole slab. In the convective 
motion the losses are localized in the layer around the URCA shell radius r^: 

P 

f* ~\~ ^conv f <i ^conv) ^conv — Ctp \j p ' (-'^"^) 



here ^conv is taken from the mean free path model. The local rate is obtained from (p^), if we 
take into account that the whole heating is concentrated inside the layer ([l9|). We get than 



= (20) 

Convective velocity suffers from additional damping due to URCA shell, because kinetic 
energy of the convection is a source of both nonequilibrium heating and of additional neutrino 
losses. The only relation of the convetional mixing length model of the convetion (see e.g. 
P]) should be modified, determining the convective velocity with an additional damping 



g-^, (21) 

P ^conv 



where Emnv = 4(5;/ is found from ([l^). These relations may be applied for description of 
URCA shell convection only for sufficiently strong convective motion, when the first term in 
( |2l| ) exceeds considerably the second one. The equation ( [2l| ) has roots only when 
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< sTp (If) 

Violation of this inequality may result in an abrupt termination of the convection in the layer 
([T9I ) around the URCA shell. When analysing the URCA shell convection in star, it would 
be premature to predict the results of evolutionary calculations with account of convective 
URCA shell before such calculations are done. Two possibilities may be expected . One 
is connected with obtaining of a definite result which has a little sensitivity to the input 
parameters of the problem, such as Op, -Fii/2, accepted rates of nuclear reactions, neutrino 
losses etc. Another possibility could be a great sensitivity of the result to the same input 
parameters. If the second possibility would be realized we could still remain in situation 
of ambiguity, because the set of the input parameters for presupernovae model cannot be 
established with a sufficient precision. 
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